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We present recent results in precision multiboson (+jet) phenomenology at the LHC. Results 
for diboson + jet, triboson, and also for W ±r yy+ jet will be discussed focusing on the impact 
of the perturbative corrections on the expected phenomenology. 

1 Introduction 

Processes with multiple electroweak bosons are important channels to test the Standard Model 
(SM) at the LHC. They are important backgrounds to SM and also to beyond standard physics 
searches. As a signal, they allow us to obtain information on triple and quartic couplings, and 
therefore, to quantify deviations from the SM prediction through, e.g., anomalous couplings. 

To match the experimental accuracy, precise and reliable predictions beyond the leading 
order (LO) perturbative expansion are required not only for cross sections but also for differential 
distributions. As part of such a program we have, in the past, determined next-to-leading- 
order (NLO) QCD corrections for the production cross sections of all combinations of three 
electroweak bosonsOH, to IFy + jetffl , WZ + jetffl and also to W ±/ yy-\- jet^, available in 
the VBFNLO package®. In all cases, leptonic decays of the electroweak bosons were included 
in the calculations. For the production of three weak bosons and also for TF ± 77, these results 
were verified against independent calculations! 9 ! 10 ^ which are available for on-shell bosons and 
neglecting Higgs boson exchange. 

In these proceedings, we review results for TF ± 77, VF ± 77 + jet and W^y/W^Z +jet, 
including their leptonic decays and full off-shell effects, in Section [2j [3] and [4j respectively, for 
the LHC at 14 TeV. We summarize in Section [5l 

2 ^77 

Among the triple vector boson production channels, TF ± 77 has turned out to be of particular 
interest. H /± 77 production is sensitive to the WWy and WWyy vertices^. In addition, a 
final state with two photons and missing transverse energy is relevant in a variety of beyond 
the standard model scenarios^ in gauge-mediated supersymmetry breaking, for instance, the 
neutralino is often the next-to-lightest supersymmetric particle and decays into a photon plus a 
gravitino, giving a signal of two photons and missing Et- In Ref.^, a study of the backgrounds 
for supersymmetry motivated di-photon production searches has been performed, pointing out 
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the relevance of the W^jj production process as a SM background in case of electron misiden- 
tification. Another possible application is an estimate of backgrounds when searching for WH 
production, followed by Higgs decay to two photons. 

We compute the NLO hadronic cross section by straightforward application of the Catani- 
Seymour dipole subtraction^. The loop contributions are evaluated using the Passarino- 
Veltman scheme up to four-point functions^ and the Denner-Dittmaier reduction^ for five 
point integrals and we perform various cross checks to validate our implementation, Refs.lHllS 
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Figure 1: Example of the three topologies contributing to pp — > ii/77 + X 

The NLO virtual corrections result from one-loop diagrams obtained by attaching a gluon 
line to the quark-antiquark line in diagrams like the ones depicted in Fig. [T] We combine the 
virtual corrections into three different groups, which include all loop diagrams derived from 
a given Born level configuration. This leaves us with three universal building blocks, namely 
factorizable corrections (Virtual-born) and corrections to two (Virtual-box) or three (Virtual- 
Pentagons) vector bosons attached to the quark line. For our numerical results, we use the 
CT10 parton distribution set 113 with a s (m z ) = 0.118 at NLO, and the CTEQ6L1 set ESI with 
<Xs(.mz) = 0.130 at LO. We impose a set of minimal cuts on leptons, photons and jets, namely, 



p Te{l) > 20 GeV |y £(7) |<2.5 R~ n > 0.4 R^ > 0.4 R je > 0.4 R j7 > 0.7 (1) 
as well as an isolation criteria a la Frixione^ for the photons, 

E< E Ti 9(5- R tl ) < p T ^ 1 - COs6 x (for all 6 < S ) , 

1 — COS Ofj 

where 5q is a fixed separation which is set to 0.7. We consider decays to the first two lepton 
generations, i.e., W — > eu e (+~f), //^(-py) and these contributions have been summed in Fig. [2j 
where we show numerical results for W + jj production within the cuts of Eqs. (|TJ [2]) . On the left 
panel, we show the overall scale variation of our numerical predictions at LO and NLO: the NLO 
K-factor is large both in absolute value (~ 3) and compared to the LO scale variation. The 
NLO scale uncertainty is about 10% when varying the factorization and the renormalization 
scale fi = fj,p = (J>r up and down by a factor 2 around the reference scale no = mvK 77 and 
is mainly driven by the dependence on fj>R- The large size of the NLO corrections partially 
originates from new gluon induced channels entering first at NLO, gq — > W ±r yj q, which are a s 
suppressed, but enhanced by the large gluon pdfs at the LHC. Since these 1-jet contributions 
to the 0(a s ) cross section are only determined at LO, and are unbalanced against the virtual 
part, their scale variation is large. In fact, most of the scale variation of the total NLO result is 
accounted for by the real emission contributions, defined here as the real emission cross section 
minus the Catani-Seymour subtraction terms plus the finite collinear terms. This is more visible 
in the right panels, where we show the scale dependence and compare the size of the different 
parts of the NLO calculation. As for the relative size of the NLO terms, the real emission 
contributions dominate and are even larger than the LO terms plus virtual terms proportional 



(2) 




Figure 2: Left: Scale dependence of the total LHC cross section for pp — » 77 + |J T + X at at LO and NLO, 
within the cuts of Eqs. [ 1|2 1. The factorization and renormalization scales are varied in the range from 0.1 ■ fio to 
10 • /j-o- Right: Same as in the left panel but for the different NLO contributions at fj,F = fiR = £/xo- 



to the Born amplitude. Non-trivial virtual contributions, namely the interference of the Born 
amplitude with virtual-box and virtual-pentagon contributions, represent less than 1% of the 
total result and their scale dependence is basically flat. In the left panels, we also show results 
for additional jet veto cuts, requiring pTj < 50 GeV or pxj < 30 GeV. While it is evident that 
the renormalization scale variation is highly reduced by a jet veto, this reduction should not be 
interpreted as a smaller uncertainty of the vetoed cross section: a similar effect in W^jj and 
WZj and W ±r yjj production could be traced to cancellations between different regions of phase 
space and, thus, the small variation is cut-dependent^ as shown in the following sections. 

Among the triple vector boson production channels, W^'yy production is the one with the 
largest K-factor for the integrated cross section. In RefP (see also 1111 ), it was shown that this is 
due to cancellations at LO driven by a radiation zero^. The radiation zero at NLO is obscured, 
similar to W ±r f production^, by additional real QCD radiation, VF ± 77 +jet, as part of the 
NLO contributions. An additional jet veto-cut might help in the detection of the radiation 
zero, while reducing also the scale uncertainties for the relevant distributions. However, this 
procedure raises the question of the reliability of the predictions due to the aforementioned 
problem with the exclusive vetoed samples. Furthermore, the remaining scale uncertainties at 
NLO QCD are due to unbalanced gluon-induced real radiation computed at LO, e.g., gq — > 
H /± 77 q. To realistically asses the uncertainties, also concerning anomalous coupling searches, 
and as an important step towards a NNLO QCD calculation of W^jj , we have calculated 
PF±77 +jet at NLO QCD. 



3 W ± 77+ jet 

This is the first calculation falling in the category of VVV + j, which includes the evaluation 
of the complex hexagon virtual amplitudes, which poses a challenge not only at the level of 
the analytical calculation, but also concerning the CPU time required to perform a full 2—7-4 
process at NLO QCD. 

For the virtual contributions we use the routines computed in Ref. At the numerical 
evaluation level, we split the virtual contributions into fermionic loops (Virtual-fermionbox) 
and bosonic contributions with one (Virtual-box), two (Virtual-pentagons) and three (Virtual- 
hexagons) electroweak vector bosons attached to the quark line. This procedure allows us 
to drastically reduce the time spent in evaluating the part containing hexagon diagrams as 



explained in Refs.HEH The numerical stability of the hexagons' contributions is discussed in 
detail in RefM 
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Figure 3: Scale variation of the ^z/77+jet production cross sections at the LHC {I = e, fi). The cuts are 
described in the text and we choose fiR = /if = mw-y, as central dynamical reference scale. The right panel 
shows the individual contributions to the NLO cross section according to our classification of topologies. We also 
show results where we have applied a veto on events with two identified jets having both p^, > 50 GeV 



We use the same input parameters as for VF ± 77 production and apply the cuts of Eqs. (Tj2). 
Further details on the parameter choices can be found in Ref.^. Again, we consider decays 
to the first two lepton generations, i.e, W — > eu e {+^), ^^(+7) and these contributions have 
been summed in Fig [3] and |4j 

We compute total K factors of 1.43 (1.48) for VF + 77+jet (W~77+jet) production at the 
LHC, values which are quite typical for multiboson+jet production as found in Refs. EMU! and 
partially originated by new e.g, gg and qq induced channels. This moderate K-factor (as com- 
pared to corrections of ~ 300% for W^'yj production) indicates, as expected, that the W^jj + 
jet production channel is not affected by radiation zero cancellations. 
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Figure 4: Differential maxp^ and mw-y-y distribution for inclusive and exclusive Z + ^77+jet production. 



The scale dependences of the W + 'yyj and W 77J production cross sections turn out to be 



modest: when comparing fJ>R = /J>f = ^mw-n f° r £ = 0.5 and £ = 2, we find differences of 10.8% 
(12.0%), respectively, see Fig. [3} 

The phase space dependence of the QCD corrections is non-trivial and sizable (we again 
choose = hf = myi/ 77 ). Vetoed real-emission distributions are plagued with large un- 
certainties (FigjiJ left) — a characteristic trait well-known from yy+jet phenomenology HEH 
Additional parton emission modifies the transverse momentum and invariant mass spectra in 
particular. The leading jet becomes slightly harder at NLO as can be inferred from the differen- 
tial K factor in the bottom panel of Fig. |4j When comparing precisely measured distributions 
in this channel against LO Monte Carlo predictions, the not-included QCD corrections could be 
misinterpreted for anomalous electroweak trilinear or quartic couplings l 4 l 6 l 25 l arising from new 
interactions beyond the SM. 



4 W ± j/W ± Z + jet 

NLO corrections to pp -^-W^^f/W^ 1 Z +jet cross section have been computed in Refs.0H, and in- 
cluding anomalous couplings in Refs.SH All off-shell effects were included. Similar observations 
as in VF ± 77 + jet were found. When varying the factorization and renormalization scale by a 
factor 2 around fixed values of hq=W0 GeV, one finds modest scale variations. Vetoed samples 
pick up large uncertainties (FigjEJ left). K- factors are around 1.4 at the LHC and they vary over 
phase space. Two examples are shown for these processes in Fig. [5j including the sensitivity to 
anomalous couplings in the pj, differential distributions for W^j +jet production for different 
choices of anomalous parameters (An, ko)- 
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Figure 5: Left: Differential distribution for pT,min,l f° r inclusive and exclusive W Z+jet production. Right: 
Sensitivity to anomalous couplings for l~ + jet in the p^, distribution. 



5 Summary 

The QCD corrections for vector boson production in the diboson + jet, triboson and triboson 
+ jet channels are large and exceed the expectations driven by LO scale uncertainties. Total K- 
factors up to 3 for W ±r y^ have been reported. The size of the QCD corrections for W 7 /W^ Z + 
jet and W^j^f + jet production is around the 40% level. Corrections can be larger for differ- 
ential distributions and therefore have to be considered for a precise comparison of data to SM 



predictions for all these processes. Finally, we have shown that the diboson + jet production 
channels are sensitive to anomalous coupling searches through differential distributions. 
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